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p_; Abstract 



Observation of the ultra-high energy neutrinos, in particular, detection of v T 
from cosmologically distant sources like active galactic nuclei (AGN) opens 
new possibilities to search for neutrino flavor conversion. We consider effects 
of violation of the equivalence principle (VEP) on propagation of these cosmic 
neutrinos. Two effects are studied: (1) the oscillations of neutrinos due to the 
VEP in the gravitational field of our Galaxy and in the intergalactic space, 
(2) the resonance flavor conversion driven by the gravitational potential of 
the AGN. We show that the ultra-high energies of the neutrinos as well as 
cosmological distances to the AGN, or strong AGN gravitational potential 
will allow one to improve the accuracy of test of the equivalence principle by 
21 orders of magnitude for massless or degenerate neutrinos (A/ ~ 10~ 41 ) 
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and at least by 12 orders of magnitude for massive neutrinos (A/ ~ 10~ 28 x 
(Am 2 /leV 2 )). Experimental signatures of the transitions induced by the VEP 
are discussed. 
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I. INTRODUCTION 



Cosmologically distant objects such as the active galactic nuclei (AGN) can be intense 
sources of high-energy neutrinos 0. The flux of these cosmic neutrinos is flavor non- 
symmetric. Models predict that the r - neutrino flux is at least 3 orders of magnitude 
smaller than the fluxes of the electron and muon neutrinos. This opens a unique possibility 
of searching for effects of neutrino flavor transitions over intergalactic distances. 

It was suggested recently [[J that large deep underwater and ice neutrino detectors will 
be able to identify the events induced by r neutrinos with energies above 1 PeV. The 
r-neutrinos produce the characteristic double-bang events. The first bang comes from the 
charged current interaction of the r neutrino and the second one originates from the hadronic 
decay of the r-lepton. At PeV energies the r-lepton tracks have typical lengths of the order 
of 100 m and the two bangs are clearly separable. Moreover, a selection criterion of greater 
energy of the second bang in comparison with the first one makes the events to be essentially 
background-free 0, thereby guaranteeing the unambiguous detection of the r-neutrinos. 

It was estimated that the experiment is sensitive to the transition probabilities greater 
than 

P > (3 - 5) • 10~ 3 . (1) 

The transition can be induced by the vacuum oscillations 0. Apart from the appearance 
of the relatively large z/ r -flux one expects also a modification of number of the v e and 
events. Cosmological distances to the AGN will allow one to probe the range of oscillation 
parameters: Am 2 > 10 _15 eV 2 , sin 2 Id > 2P « 10~ 2 . In particular, strong effect is expected 
for the values of parameters needed to solve the atmospheric neutrino problem ||. 

In this paper we will consider the flavor transitions of the high-energy cosmic neutrinos 
induced by possible non-universality in the gravitational couplings of neutrinos. We show 
that the ultra-high neutrino energies and cosmological distances, or strong AGN gravita- 
tional potential, conspire to lead to great sensitivity in exploring the tiny non-universality. 
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The non-universality of the gravitational couplings of neutrinos can lead to the neutrino 
flavor oscillations |3Q. Search for the neutrino oscillations in accelerator neutrino experiments 
and solar neutrino observation turns out to be a powerful tool of exploring possible violation 
of Einstein's equivalence principle PHT2f . Moreover, the oscillations and neutrino conversion 
induced by violation of the equivalence principle supply viable mechanism which could solve 
the solar and the atmospheric neutrino problems simultaneously • 



The sensitivities to non-universality achievable in these methods vary depending upon the 
experimental means, on the presence or absence of matter effect, and on values of neutrino 
masses. Existing neutrino accelerator experiments restrict the non-universality parameter 
Af up to the level of 1CT 15 [0,|TT],|12| , assuming that the supercluster's gravitational potential 



is of order 10 -5 (see below.) It will be possible to improve it in the planned long baseline 
accelerator experiments which can reach the sensitivity 10~ 18 ||,[7|. 

In the case of massless (or degenerate) neutrinos studies of the solar neutrinos can be 
sensitive to 10 -20 fll^f . This limit corresponds to the vacuum long length oscillations. In the 
case of the resonant flavor conversion of massive neutrinos due to the matter effect in the 
solar interior the sensitivities are at most 10~ 15 — 10 -16 P JiTl|JT^] . 

It has also been discussed that for Af > 4- 10~ 14 the conversion induced by VEP leads to 
substantial effect on supernova dynamics [fCfl . The arrival time difference between neutrinos 
and photons from SN1987a gives a modest bound of the order of 10~ 3 |I3| . 

We may conclude that by using any neutrino sources and any experimental devices so 
far considered it appears difficult to go far beyond the sensitivity of ~ 10~ 20 . 

In this paper we will show that observation of the ultra-high energy neutrinos from 
cosmologically distant sources will allow one drastically improve the accuracy of testing the 
equivalence principle. The paper is organized as follows. In Sec. II we will consider the 
gravity- induced oscillations of massless (or degenerate) neutrinos and estimate a sensitivity 
to violation of the equivalence principle (VEP). In Sec. Ill the gravity effects in the presence 
of nonzero neutrino masses and vacuum mixing are discussed. In particular, we describe the 
resonance flavor conversion driven by the gravitational potential. In Sec. IV we will apply 



the results to neutrinos from the AGN. In Sec. V the experimental signatures of the VEP 
effects are considered. In Sec. VI we summarize the results. 



II. GRAVITATION ALLY INDUCED OSCILLATIONS OF NEUTRINOS FROM 

THE AGN 

Let us restrict ourselves to the two-flavor case, (z/ M , z/ r ), for simplicity. According to the 
hypothesis of violation of the equivalence principle (VEP) Q the flavor eigenstates and 
v T are the mixtures of the gravity eigenstates, u 2g , z/ 39 , whose gravitational couplings f 2 G 
and f 3 G, where G is the Newton constant, are different f 2 7^ f 3 - Evidently, fa ^ 1 at least 
for one neutrino. Introducing gravitational mixing angle, 9 g , one can write 

Vy. = cos O g v 2 g + sin 6 g v 3g , v T = - sin Q g v 2g + cos 9 g u 3g . (2) 

The non-universality of the gravitational couplings can be parametrized as 

The v 2g and v 3g neutrinos feel gravitational fields with slightly different strengths. This 
leads to a difference in the energies of the eigenstates (i.e. to the level energy splitting): 

V g = ±AfE*(x), (4) 

where E is the energy of neutrino, and $(x) = MG/r is the gravitational potential at 
distance r from an object of mass M in the Keplerian approximation. The energy level 
splitting induces a relative phase difference between wave functions of v 2g and u 3g which 
results in neutrino flavor oscillations in the same fashion as in the mass-induced case. 

Let us first suggest that neutrinos are massless or have equal masses. As we will show 
later the matter effect is negligibly small. In this case the propagation of neutrinos has a 
character of oscillations with the depth fixed by 9 g and with the length, l g , determined by 
V 9 : 
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The oscillation probability can be written as 

P = sin 2 2fl 9 sin 2 (a/E J $(x)dx^j . (6) 

The distinctive feature of (f|) is that the oscillation length is inversely proportional to the 
neutrino energy, in contrast with the vacuum oscillation length which is directly propor- 
tional to E. Therefore the sensitivity to Af increases with the neutrino energy and with 
the path-length in the gravitational field. Both factors are present for neutrinos from the 
AGN. The AGN are believed to produce high energy neutrinos with spectrum extended to 
E ~ 10 PeV, m and typical distances from our Galaxy to the AGN are L^gn ~ 100 Mpc. 



Let us estimate a sensitivity to Af. For this purpose we should find the integral 
I = J Q(x)dx along the neutrino trajectory which appears in (||). The intergal has three 
contributions: 

/ = I AGN + I IG + I G , h = J ^(x)dx, (i = AGN, IG, G) , (7) 

where &(x)agn, §(%)ig, an d §(x)g are the potentials created by the AGN itself, by all 
bodies in the intergalactic space and by our Galaxy, respectively. 

We get Iagn^) ~ — l-Rslog {r/R e ) for a radial trajectory, where Rs is the Schwarzschild 
radius of AGN: R s ~ 3 x 10 n (M AGA r/10 8 M o ) m and R e w (10 - lO 2 )^ is the radius of 
the neutrino emission region. Using M^gn ~ 10 8 M Q as a typical mass of the AGN, we find 
Rs ~3 x 10 11 m and R e ~ 10 13 m , and consequently 

Iagn ~ 10- 10 Mpc. (8) 

In the intergalactic space the effect is dominated by the gravitational field of the so 



called Great Attractor [Q. This supercluster is located at the distance (43.5 ± 3.5) 
Iiq 1 Mpc from the Earth, where ho is in the range 0.5-1.0. The mass of the Great At- 



tractor is about M sc ~ 3 x 1O 16 VM , where M Q is the solar mass. |15| |. In the Kep- 
lerian approximation the gravitational potential of this supercluster can be estimated as 
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= -5.2 x lO- 6 (i2/lOOMpc) _1 (M sc /lO 16 M ). (Note that the weak-field approximation 
still applies). Then the integral Ijq satisfies the inequality 

for any trajectory of length L within radius of 100 Mpc. 

The contribution of our Galaxy equals to Iq = —GMq log (Lagn /t) for radial trajectory 
Using Mq ~ lO n M and r = 10 kpc for the mass and the radius of the Galaxy, respectively, 
and L AGN = 100 Mpc we get I G ~ -10~ 7 Mpc. 

The supercluster gravitational effect dominates: I ~ Ii G . The reason is that the path- 
length in the gravitational field tends to cancel the inverse distance dependence of the 
Keplerian potential. Therefore, the oscillation probability is essentially governed by the 
mass of source of the gravitational field and the mass relation M sc 3> M G ^> Magn leads 
to I sc > I G > Iagn- 

It is straightforward to make an order-of-magnitude estimation of the sensitivity. For 
E— 1 PeV and I IG (at $ = 10~ 5 , and the distance L = 100 Mpc) one gets 

\h G E\ > \$EL\ = 1.5 x 10 41 (10) 

which means according to (^) that the phase of oscillation of order unity will be obtained 
for Af > 10" 41 . 

If the great attractor is a fake object, then the dominant effect would be due to the 
gravitational field of our Galaxy: I w Iq. In this case the expected sensitivity to Af is 
~ 10 -37 which still implies an improvement by more than 17 orders of magnitude. This Af 
can be considered as the conservative estimation of the sensitivity. 

Let us show that in spite of cosmological distances the matter effect on the neutrino 
conversion can be neglected. Consider the u e — v T system for which matter effect appears 
in the first order in the weak interactions. The cosmological baryon density estimated from 
the nucleosynthesis is p% ~ 10~ 31 g/cm 3 . This gives the width of matter in the inter- 
galactic space: die = Pb • Lagn ~ 3 • 10~ 5 g/cm 2 . According to spheroid-dark corona 
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models [OJ the matter density of our Galaxy is p ~ (1 — 10) x 10~ 25 g/cm 3 . This leads 
to the width dc ~ 3 ■ (1 — 10) x 10~ 3 g/cm 2 . Finally, the width of matter crossed by 
neutrinos in the AGN is estimated as (Iagn ~ (10~ 2 — 10 _1 ) g/cm 2 PT7[] . Thus total 
width, dtotai ~ &AGN ~ (10" 2 — lO^ 1 ) g/cm 2 , is much smaller than the effective width 
do = \f2-Km n I Gf ~ 2 • 10 9 g/cm 2 needed for appreciable matter effect. For — z/ T channel 
the matter effect appears in high order of perturbation theory and the required effective 
width is even larger. 



III. GRAVITATION ALLY INDUCED TRANSITIONS IN THE PRESENCE OF 

NEUTRINO MASSES 

Most probably neutrinos are massive and mixed. The gravitational effects themselves can 
generate via the nonrenormalizable interactions the neutrino masses of the order v 2 /Mp ~ 
10 -5 eV where v is the electroweak scale and Mp is the Planck mass. Moreover, 
there are some hints from solar, atmospheric, cosmological as well as accelerator data that 
neutrino masses are even larger than that value. Forthcoming experiments will be able 
to check the hints. In this connection we will consider the gravitational effects in the 
presence of neutrino masses and mixing, assuming that the latter will be determined from 
the forthcoming experiments. 

In the presence of the vacuum and gravity mixing the effective Hamiltonian of neutrino 
system in the flavor basis v = {y^, z/ T ) reads ||] (for notation we follow ||): 

f-C S\ f-Cg Sg\ 

H(x)=si \+V g i , (11) 

V s c J \ s g c g I 

where 5 = Am 2 /4£ with Am 2 = m^ 2 — m 2 2 , c = cos 29, and c g = cos 29 g , etc.. As we 

have shown in Sec. II the matter effects can be neglected. For antineutrinos the VEP may 

differ from that for neutrinos: Af ^ Af, so that one may expect a difference of the effects 

in neutrino and antineutrino channels. We assume at the moment that all the parameters 



S 



of the Hamiltonian ([11]) are real. In general, the complex phase can be introduced in (|TT]) 
which has some physical consequences fll2| . We will discuss effects of the phase at the end 
of this section. 



According to ( |TT] ) the mixing angle in medium, 9 m , is fixed by 

S5 + VgSg 



tan2fl m = — . (12) 

CO + VgCg 

Evidently, for V g ^> 5 the mixing is determined by gravity: 9 m ~ 9 g . And for V g <C 8 the 
mixing angle equals the vacuum angle: 9 rn « 9. 
As follows from (|I2]) the mixing angle is zero at 

K° = -5- • (13) 



s 



If 9 g = 9 one has from ( |T2"P tan26 l m = tan 29. In this case the evolution of neutrino state 
is reduced to oscillations with the constant depth, sin 2 29, and the oscillation length 

The phases due to mass difference and the gravitational effect add up. The gravitational 
phase dominates if V g > 8, or explicitly, if Af > Am 2 /(2E 2 &). For supercluster potential, 
$ sc ~ 10" 5 , and E — 1 PeV we get from this inequality Af > 5 x 1(T 36 , 5 x 1(T 28 , 5 x 10~ 25 
for Am 2 = 10 -10 , 10~ 2 , 10 eV 2 , respectively. However, in all these cases the oscillations are 
averaged and even if Am 2 will be known it is impossible to identify the gravity effects. 



If 9 7^ 9 g , the Hamiltonian (|TT]) can lead to the resonance flavor conversion due to change 



of V g with distance, i.e. to the gravity version of the Mikheyev-Smirnov-Wolfenstein (MSW) 



mechanism ||19|| . Note that apart from a brief remark in the consideration in the existing 
literature is restricted to the level crossing driven by matter density change with distance. 
Here we describe the level crossing driven by gravitational potential change rather than 
matter density change. 



As follows from (|Tl|) the resonance condition (condition for maximal mixing) is 

V. = ~ , (15) 

and the resonance value of the potential equals 

Am 2 cos 2^ 

* = ~2E 2 Af cos20 s ' ( } 

Correspondingly, Vg in the resonance is V g R = AfE^ R /2. 
The adiabaticity condition in the resonance reads |fl,[T0|| 

(5sin20 + V g sin20 9 ) 2 



f^cos2^ 



> 1. (17) 



resonance 



This condition simplifies under the Keplerian approximation. Substituting in (]17|) dV g /dr 
—Vg/r and using the resonance condition (|15D we can rewrite fll7|) as 



(tan2#-tan2#„) 2 cos2#„ > — , (18) 

r R 5 

where r R is the radius at which the resonance condition is fulfilled: r R = r R (Af). 

The conditions (15, [T8|) determine the sensitivity regions of parameters Af, 9 g . Indeed, 



for fixed Am 2 and c ~ c g ~ 1 the minimal and maximal values of the potential $ determine 
via the resonance condition the range of Af for which the resonance neutrino conversion 
may take place. Maximal value of r R at which the resonance condition is fulfilled gives the 
lower bound on mixing angles through the adiabaticity condition flUf). If the adiabaticity 
condition is satisfied with vacuum mixing angle alone, i.e., sin 2 29 « tan 2 29 3> (r^) -1 , 
then 9 g can be arbitrarily small. In this case the role of the gravitational effect is just to 
split levels. On the contrary, for 9 = one has the lower bound on the gravitational mixing 

sin 2 2^>^-. (19) 
r R 5 

For fixed mixing angles the adiabaticity condition can be rewritten as the lower bound 



on Af. Indeed, substituting 5 in the resonance condition ( flq) into (|18D we find 

A f » 2 >< 10 ' 33 ^ (M*™) ~ l (20 ) 

J (tan2#-tan2fl 9 ) 2 cos2# 9 UPeV/ V lo8M ©/ ' 
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If the adiabaticity condition is fulfilled, then the transition probability is determined by 
the initial and final values of mixing angle: 

P a = ~ (1 - cos 29 mi cos 26 mf ) . (21) 

In this connection let us consider a dependence of mixing angle 9 m on the potential <3> or 
level splitting V g = AfQE/2. The angle 9 m as the function of $ crucially depends on the 
sign of Am 2 Af, and on whether 9 g > 9 or 9 g < 9. We focus first on the resonant channel, 
Am 2 Af > 0. 

1) . 9 g < 9. In this case one has \V g °\ > \V g R \, i.e. the value of potential which corresponds 
to zero mixing fll3|) is bigger than the resonance value. For the initial splitting \V g \ \V g R \, 
the gravity mixing dominates and 9 m ~ 9 g + tt/2. With diminishing \V g \, the angle 9 m 
decreases and mixing becomes zero (9 m = tt/2) at V g = V g ] then 9 m crosses the resonance 
value, 9 m = ir/4, and for \V g \ <C \S\ approaches 9 m = 9. 

2) . 9 g > 9. Now \V g °\ < \V g R \. At |Vj| > \V g R \ one has 9 m ^9 g - tt/2. With diminishing 
\V g \ the angle 9 m increases and crosses resonance value 9 m = — ir/4. At V g = V g ° the angle 
9 m vanishes so that sin 2 29 m = 0, and then 9 m approaches the vacuum value 9. 

In the nonresonant channel, movement of the angle 9 m is simpler. Under the same vari- 
ation of V g as above it starts from 9 m = 9 g and ends up with 9 without crossing the points 
of zero mixing and resonance irrespective of the relative magnitudes of 9 g and 9. 



Suppose that the initial potential (the potential at the production point) is much larger 
than the one at resonance, $j 3> $_r, and the final potential is much smaller than the value at 
resonance, -C In this case the transition probability in the adiabatic approximation 



(pID becomes: 



P a = ^(I±cos20 fl cos20), (22) 
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where the plus sign is for the resonance channels (Am 2 Af > 0), and the minus sign is for 
the non- resonant channels (Am 2 Af < 0). 

Let us mark one interesting feature related to zero mixing at V g — V g ° (|1"3|). If the initial 
(final) potential is such that V g = V g ° for 6 g < 9 (6 g > 9), then the transition probability in 
the resonant channel reduces to P = cos 2 9 (P = cos 2 9 g ) as in the case of flavor conversion 
in matter. 



Let us discuss finally effects of possible complex phases in the Hamiltonian (|TT1). After 
redefinition of the neutrino wave function only one phase survive which can be put in the 
gravitational term of the Hamiltonian flTT|) [ |i~2]j : 



V 



c g s g e 



-i2a 



s g e 



r2a 



(23) 



The angle a is the relative phase of the vacuum and gravitational contributions to the 
mixing (nondiagonal elements), and it is present when both of these contributions exist. Let 
us consider physical effects of the phase. 

The Hamiltonian (|ll|) can be rewritten as 

-cS - V g c g Pa e- i2 ^ 



H(x) 



Pa e i2 ^ c5 + V g c g 



where 



Pa 



with £ = 5s/V g Sg, and 



5s + V g s g e- i2a 



VgSg V e + 2£cos2a + l 



tan 2^- ^ Sin2a 



5s + V g s g cos lot 



Additional redefinition of the fields 



(24) 



(25) 



(26) 



^ = e*%, i/ T = e-*>v T , 



(27) 



leads to the evolution equation for (y' u' T ) with the Hamiltonian 
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-c8-VgCg-lf> Pa 

H(x) =| . 1 , (28) 

Pa C5 + VgCg + 1p , 



where ip = dif) /dr. 



From this we get the following consequences. 



1) . There is no complex phases in the Hamiltonian fl28|). Thus, the presence of the phase 
a does not lead to the CP— or T— violating effects in the two generation it was 
expected from the beginning. 

2) . The resonance condition is modified: 

c8 + V g c g + i> = . (29) 

The shift of the resonance, if), is absent in the case of constant gravitational potential for 
which if) = 0. If the adiabaticity condition is fulfilled in the absence of a then the shift of 
the resonance position is negligibly small. Indeed, for a maximal phase sin 2 2a = 1, we find 
from ( pSD that ip < ji- and then after suitable modification including the effect of phase, 
the adiabaticity condition leads to if) <C V g s g - The latter condition means that the shift of 
the resonance is much smaller than the width of the resonance. 

3) . Nonzero phase a modifies the dependence of the mixing parameter sin 2 29 m on V g . The 
mixing is proportional to the nondiagonal element of the Hamiltonian (p8|) p a . For £ ^> 1 
or £ 1 it follows from ( p5|) that p Q ~ p , where p corresponds to the case of zero a. That 
is, the effect of a is small when the gravity or vacuum effect dominate. 

The strongest influence of a is for £ = — 1 which corresponds to the point of zero mixing: 
V g = Vq. In fact, the phase remove zero mixing at Vo- 

4) . If £ cos 2$ > 0, then p a < p ; the phase leads to decrease of mixing and to narrower 
resonance peak. For £ cos 28 < we get p a > po, i.e., the phase enhances mixing. 

5) . The presence of phase affects also the adiabaticity condition. However, it can be shown 
that if the adiabaticity is fulfilled at a = 0, this influence is weak apart from some exceptional 
values of the phase. 
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Thus we see that in many physically interesting situation the role of the complex phase 
is quite small even if a ~ 0(1), and in what follows we will discuss the case a = 0. 



IV. FLAVOR TRANSITION OF THE AGN NEUTRINOS 

Let us consider the flavor transitions of neutrinos from the AGN using the results of Sec. 
III. Following the scenarios described in Q, we assume that neutrinos are produced within 
the region located at the distance R e = (10 — 100)Rs from the center of AGN, where R$ 
is the Schwarzschild radius: R$ — 3 x 10 u (Magn/^ 8 M q ) m. For radii larger than the 
neutrino production point we may use the Keplerian approximation for the potential of the 
AGN. The total potential probed by neutrinos on the way to the Earth is 

$(0 « $4cv(r) + $/ G = $agn (^r) + ®IG ■ (30) 

Here, 

is the AGN potential at the neutrino production point and, for simplicity, we take the po- 
tential in the intergalactic space to be constant: $jg — 10~ 5 . Therefore, at the neutrino 
production point the potential §agn dominates over the supercluster and the galactic po- 
tentials. For AGN located at 100 Mpc from us it is about 3 and 7 orders of magnitude larger 
than the potentials of the great attractor and the Milky Way Galaxy, respectively. In what 
follows we neglect the potential of our Galaxy. 

For fixed A/ the dependence of the transition probability on the neutrino energy (or 
E 2 /Am 2 ) is the following. For small energies the mass-induced vacuum oscillation effect 
dominates and 

1 F 2 1 

Pas- sin 2 26 for — — < — = — . (32) 

2 Am 2 2$ AGN Af 1 ; 
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For larger energies the transition probability equals P a in (0), if the adiabaticity condition 
is fulfilled: 

p ~ p - for ^>M^A7' (33) 

Moreover, in the region where AGN potential dominates over intergalactic potential the 

resonance conversion takes place. This corresponds to 

1 E 2 1 

2$° AGN Af < Am 2 < 2$ /G A/ " ( } 

For these energies the transition probability is larger than 1/2, and can be close to 1. 

The weaker the potential of the supercluster the wider the resonance-effective region of 

parameters. Moreover, if the $/g < 10 -7 the resonance conversion may take place in the 

gravitational field of our Galaxy too. 

For higher energies the mass splitting can be neglected and the dominant effect is the 

one due to gravitationally induced oscillations: 

P « - sin 2 26, for — — > — . (35) 

2 9 Am 2 2§ IG Af v ; 



Once the adiabaticity condition is satisfied in the AGN gravitational field the transition 
probability (|22|) applies also to the nonresonant channels. 

The resonance flavor conversion is effective thanks to the fact that the mass of the AGN 
is concentrated mainly in a small region of space (~ 10~ 5 kpc) of the central black hole. 
Since our Galaxy is more massive than the typical AGN one might expect that neutrinos 
converted at the AGN could be reconverted by our Galaxy's gravitational field. It does 
not occur for ~ 10~ 5 : the gravitational potential of the Galaxy is at most $ G ~ 10~ 7 
because its mass is distributed over the ~ 10 kpc region. 

Let us find the sensitivity to Af assuming that future underwater /ice experiments will 
be sensitive to the flavor transition with probability as small as in ([[]) [Q. We can distinguish 
three ranges of Af corresponding to three energy regions defined in (^), (0) and fl3~5|). 
Region I: 
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Af < Af AGN , (36) 

where 

Here E is the average energy of the detected neutrinos. The gravity does not play any role 
and the effect of mass-induced vacuum oscillation dominates. 
Region II: 

AJagn < Af < Af IG , (38) 

where 

In this region neutrinos undergo resonance conversion driven by the potential of the AGN, 
and the transition probability can be close to 1, if both 9 and 9 g are small. 

Suppose that the factor (tan 29 — tan 29 g ) 2 cos 29 g in the adiabaticity condition ( pH|) is of 
order unity. Then the adiabaticity condition is satisfied for Af > 2x 10~ 33 . Comparing this 
inequality with fl39|) we find that for Am 2 > 10~ 4 eV 2 the adiabaticity holds in whole region 
where the resonance condition is met. If 10 -7 eV 2 < Am 2 < 10~ 4 eV 2 the adiabaticity region 
partially overlaps with the resonance region. If Am 2 < 10 -7 eV 2 there is no overlapping 
between the two regions and the resonance conversion is not efficient. 
Region III. 

A/ » Af IG . (40) 

Here gravitational effect dominates and neutrinos oscillate due to the VEP. The probability 
converges to P = | sin 2 29 g . 

V. OBSERVATIONAL SIGNATURE OF GRAVITY EFFECTS 



The underwater /ice detectors will be able to measure the ratio of number of events 
induced by v T and u^, v T jv^ as well as the ratio of v e and events, v^jv^. 
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If neutrinos are massless or degenerate, the excess of the r neutrino signal in underwa- 
ter installations will testify for the neutrino flavor transition induced by the VEP. As we 
discussed in sect. 2 the appropriate range of parameters is Af > 10~ 41 and sin 2 29 g > 0.01. 

How can one prove this? 

Evidently, the anisotropy of Vrjv^ and z/ e /V M ratios correlated to the position of the 
supercluster and/or the center of our Galaxy could be a signature of the gravitational origin 
of the excess. Let us compare the gravitational effect in the direction towards and away from 
the supercluster. Consider the AGN's at the distance R from the center of supercluster in 
the cone directed from the Earth toward and outward the supercluster. Neutrinos from such 
AGN will acquire, on the way to the Earth, the gravitational phases 

^toward = + , <p a ™ a y = , (41) 

where 



A0 = AfEM sc G 



21n- + l 

r 



, = AfEM sc G 



R 
111 1. 



(42) 



Here d is the distance from the Earth to the supercluster and r is the size of the supercluster. 
The difference of phases is described by A<p. The effect is determined by integration over R. 
(If we neglect the expansion of the Universe the flux of neutrinos from different R will be 
the same in the case of uniform distribution of the AGN). The anisotropy will be observable 
provided that the following conditions are met: 

1) . Acj) ~ 1, i.e. the phase should be large enough. 

2) . (f)(R = ct u ) < A0, where t u is the age of the Universe. This condition ensures the 
absence of the spatial averaging. 

3) . Events within sufficiently small energy interval AE < E are selected to avoid the 
averaging over the energy. 

For d — 50 Mpc, r = 1 Mpc we get A0/0 ~ 2. It means that for the supercluster the 
conditions 1). and 2). can be fulfilled for certain values of E and Af. However, observation 
of the asymmetry will require high statistics, since one should select the events within small 
solid angle, ~ (r /d) 2 , and an energy-cut is required. 
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The conditions 1). and 2). are not satisfied for our Galaxy as the source of the asymmetry. 

If an excess of v T events will not be found, then this will allow one to exclude the region 
of parameters, sin 2 29 g > 1CT 2 and Af > Af^GN, where in the expression ( pTf ) for AJagn 
one should take as Am 2 the upper experimental bound on neutrino mass difference. 

As we have discussed in Sec. Ill most probably neutrinos are massive and mixed, and 
moreover, there is a good chance that forthcoming experiments will measure Am 2 and 
sin 2 29. In the case of massive neutrinos a signature of the gravity effect is quite different 
from that of the massless case. One should look for the deviation of the observed ratio 
^t/*V s ( as we U as v el v p) from that stipulated by vacuum oscillations. The gravity induced 
mixing can both suppress and enhance z/ T -signal. 

If Af is in Region II fl38|) the gravitational MSW effect occurs in the resonance channel. 
The resonance conversion of to v T gives rise to a large transition probability P > 1/2. 
Therefore, the observation of the ratio z/ T /z/ M 's larger than unity would provide an evidence 
for the gravitational MSW effect. 

In the nonresonant channel the gravity effect is described by the transition probability 
(f2"2|) with minus sign. For 9 g > 9 {9 g < 9), the gravity effect enhances (suppresses) the 
transition. The modification is maximal if the vacuum angle 9 is small and the gravity angle 
is large, 9 g ~ 7r/4. 

If Af falls into Region III ([|(]) the signal would mimic the one due to vacuum flavor oscil- 
lation. But, since we assume that the masses and vacuum mixing angles will be determined 
by future experiments, the difference between 9 g and 9 should show up in the measured ratio 
of v T to v^. 

Let us discuss experimental signatures and estimate the sensitivity regions for the VEP 
parameters for two probable scenarios of neutrino masses and mixing. They are suggested 
by some experimental results at hand and will be checked by forthcoming experiments. 
1. Suppose the heaviest neutrino (which practically coincides with v T ) has the mass in the 
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cosmologically interesting region: = (3 — 7) eV, so that Am 2 = (10 — 50) eV 2 in the 
case of mass hierarchy. For these values of mass the present experimental bound on mixing 
angle is already sin 2 29 < 5 ■ 10~ 3 , i.e. below the sensitivity limit (|1|). We will assume also 
that mixing with electron neutrinos is negligibly small. 

In this case the effect of mass-induced vacuum oscillation can be neglected. The result for 
the ratio v T j > 0.01 in the underwater /ice detectors would then be an indicator of the VEP 
mechanism of flavor conversion. Observation of large value of the ratio: ^ T /^> 1, would 
provide a clear signature for the gravitational MSW effect. From the resonance condition 



we find, assuming the typical values for the parameters in (|37|) and (p9|), the required region 
for Af: 

Af = (10- - 5 x 10-) (^) . (43) 

Then, for fixed Af, the adiabaticity condition gives the bound on mixing angles 

(tan2fl - tan2g g ) 2 cos2fl g » 5 x 10' 9 ( - -ALg ) (J£^ ■ (44) 

If lepton mixing is similar to the quark mixing, then one expects sin 2 29 > 10~ 3 for u T — u^. In 
this case the adiabaticity condition is satisfied by vacuum mixing alone and the gravitational 
angle can be arbitrarily small. 

Non-observation of such a signal will allow one to exclude the whole region of parameters 
in (|43|) and 9 g being not too close to 9. 

Observation of a moderately-large ratio, 0.01 < Vrjv^ < 1 would be the signal for 
one of the following three possibilities; the nonadiabatic gravitational MSW mechanism 
(in resonant channel), the nonresonant adiabatic conversion with transition probability 
P a ~ sin 2 9 g (in nonresonant channel), or oscillations in the intergalactic gravitational field 
with sin 2 29 g > 0.01. In the first case, the parameter Af should take the value around ( f43|) 
and 9 g w 9 to violate the adiabaticity which is well satisfied by vacuum mixing alone. In 
the other two cases, the region sin 2 29 g > 0.01 will be probed. 
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2. Suppose that the heaviest neutrino has a mass ~ 0.1 eV and — v T mixing is large: 
sin 2 29 ~ 0.5 — 1, so that — v T oscillations solve the atmospheric neutrino problem. As the 
result of the vacuum oscillations one predicts the ratio v T jv iJ , < 1 for neutrinos from AGN. 
The observation of larger ratio, v T /v^>l, will be a signature for the gravitational MSW 
effect. In this case the VEP parameter should be in the interval 

Af = (0.3 - 1) x 10- 30 - 5 x 10~ 28 , (45) 



assuming again the typical values of parameters in the parentheses in ( p7|) and ([39]). The 
adiabaticity condition is satisfied by large vacuum mixing angle, and the angle 9 g can be 
arbitrarily small. 

In the nonresonant channel the gravitational effect manifests in a different manner. If 
Af is in Region II (0), i.e., Af > 10~ 30 , and the gravitational mixing angle 8 g is smaller 
than the vacuum angle 9, the gravitational effect suppresses neutrino transition. Namely, 
the v T signal will be smaller than the one expected for mass-induced vacuum oscillations. 
Using ( |2"2"| ) with ^ 9 < 1 we find transition probability P = sin 2 9 instead of P = ~ sin 2 29 



in the absence of the VEP. In such a way the transition probability is reduced by 40% and 
30% for sin 2 29 = 0.6 and 0.8, respectively. 

In Region III ( flOf) the averaged transition probability equals P — \ sin 2 29 g and for small 
9 g suppression can be much stronger, so that v T signal will not be observable at all. 



VI. CONCLUSION 

1. Cosmological distances ( ~ 100 Mpc) and ultrahigh energies (~ 1 PeV) of neutrinos from 
the AGN open unique possibility to improve an accuracy of testing the equivalence principle 
by 12 - 21 orders of magnitude. 

2. For massless neutrinos VEP can induce the oscillations — v T of cosmic neutrinos which 
may lead to observable v T signals in the large underwater /ice installations. The sensitivity 
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to the parameters of the VEP can be estimated as Af > 1CT 41 and sin 2 29 g > 2 • 1CT 2 . In the 
case of the nonaveraged oscillations (Af at the lower bound) one can expect an anisotropy 
of the v T signal correlated to the position of the Great Attractor. 

3. In the case of massive neutrinos the gravitational effects due to VEP can modify the result 
of the vacuum oscillations. For certain values of the parameters neutrinos may undergo the 
resonance flavor conversion driven by the gravitational potential of AGN (or of our Galaxy, if 
the intergalactic potential is sufficiently weak). The gravitational effects become important 
if Af>Af AGN ~ 10^ 28 (Am 2 /l eV 2 ). For Af ~ (1 - 10 3 )Af AGN one may expect almost 
complete transition of to v T due to the resonance conversion. A strong observable effect 
may exist for arbitrarily small sin 2 26^. For Af 3> (1 — 10 3 )Af AGN neutrinos undergo the 
gravity induced oscillations and vacuum mixing effect can be neglected. 

4. The VEP effects can be identified if the gravitational mixing differs appreciably from 
vacuum mixing. In general, the VEP effect will manifest itself as deviation of the observed 
ratios v T j and v^j u e from those stipulated by vacuum oscillations. Of course, the latter can 
be predictable only if neutrino masses and mixing are determined by forthcoming neutrino 
experiments. For small vacuum mixing the VEP can enhance the — v T transition and the 
z/ r -signal. On the contrary, for large vacuum mixing the VEP can lead to suppression of the 
z^-signal. The ratio z/r/V M > 1 is the clear signature of the resonance flavor conversion. 

5. If deviation from vacuum oscillation effects will not be found, then one will be able to 
exclude very large new region of the VEP - parameters. 
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